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ABSTRACT 


In  recent  years,  a  large  increase  in  the  number  of  computer 
controlled  systems  and  processes  has  dictated  the  use  of  devices  with 
inherently  long  time  delays.    The  result  has  been  unstable  control 
systems.     This  thesis  develops  two  methods  of  feedforward  compensation 
which  cope  with  this  problem.     The  Pade'  method  of  polynomial 
approximation  is  used  to  formulate  a  number  of  transfer  functions  which 
can  be  used  in  the  two  feedforward  compensation  schemes.    To  choose 
the  best  transfer  function    for  a  particular  problem,  a  minimum  integral 
square  error  analysis  technique  is  developed  for  use  in  the  frequency 
domain.    Computer  simulations  are  used  to  verify  the  effectiveness  of 
the  two  compensation  schemes. 
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I.     INTRODUCTION 

In  any  control  system,  with  lumped  linear  elements,  the  output  is 
theoretically  characterized  by  a  response  which  starts  immediately  after 
the  application  of  a  transient  input.     It  is  well  known  that  this  is  not 
really  true  in  many  cases,  since  a  number  of  devices  may  have  a 
transportation  lag  or  other  delay  associated  with  them.    If  the  delay  is 
small,  it  can  usually  be. ignored.     The  system  is  then  described  by  a  set 
of  linear  differential  equations  and  the  classical  methods  of  analysis 
can  be  used  to  study  the  behavior  of  the  system.    In  recent  years,  the 
large  increase  in  the  number  of  computer  controlled  systems  and 
processes  has  dictated  the  use  of  devices  with  inherently  long  time 
delays.    Examples  of  such  devices  would  be  a  digital  computer  controlled 
system,  where  the  computation  time  could  be  several  seconds,  or  an 
automated  superfractionator,  which  could  have  a  time  delay  as  long  as 
several  hours.    A  typical  system,  with  unity  feedback,  is  shown  in  Fig.   1. 

When  the  effects  of  time  delay  cannot  be  ignored,  the  primary 
problem  is  one  of  system  instability.    This  is  easily  seen  by  looking  at 
the  Bode  diagram  shown  in  Fig.   2.     This  diagram  represents  the  open 
loop  transfer  function  for  the  system  in  Fig.   1  with 

Q   (s)  =      «™L__-  (!) 
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and 
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FIGURE   2. 
Bode   Diagram   for   Equation   1. 
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The  time  delay,  T,  was  set  equal  to  1  second.    The  upper  phase  angle 
curve  represents  the  system  without  the  time  delay.    It  can  be  seen 
that  the  effect  of  the  delay  is  to  increase  the  phase  lag  as  the  frequency 
increases.    The  amount  of  the  increase  is  equal  to  the  product  of  the 
frequency  and  the  time  delay.    In  this  particular  example,  a  very  stable 
phase  margin  of  58.5°  has  been  reduced  to  24°.    If  the  time  delay  is 
increased  to  2  seconds,  the  system  would  be  unstable.    It  can  be  easily 
concluded  that  even  a  modest  delay  can  present  formidable  stability 
problems. 

This  paper  presents  a  method  of  compensating  for  time  delay  in 
closed  loop  systems,  such  as  the  one  shown  in  Fig.   1.    The  technique 
involves  feeding  a  conditioned  signal  forward,  around  the  element 
associated  with  the  time  delay.    The  feedforward  path  may  require  either 
active  or  passive  elements  to  condition  the  signal.    The  circuit  require- 
ments will  be  developed.    This  method  is  an  offspring  of  one  referred  to 

1,2 

by  P.  S.   Buckley         as  the  Smith  Linear  Predictor.     However,  Buckley's 

use  of  the  predictor  was  not  to  stabilize  a  closed  loop  system.     He  was 
concerned  with  the  effects  of  an  external  disturbance  (noise)  being 
introduced  into  the  loop  between  the  controller  and  the  plant.     He 
assumed  the  system  to  be  stable  before  the  disturbance  was  applied.    The 
predictor  was  only  a  method  for  reducing  the  effect  of  the  disturbance. 
The  Smith  Linear  Predictor,  as  presented  in  Buckley's  paper,  is  shown 
in  Fig.  3. 

It  can  be  observed  that  the  circuit  requirements  for  the  signal 
conditioner  in  the  feedforward  path  cannot  be  synthesized  with  lumped, 
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linear,  passive  network  elements.    A  suboptimal  network  will  therefore 
have  to  be  used  in  the  feedforward  path.    The  Pade'  method  of  polynomial 
approximation  will  be  derived  and  used  to  find  an  acceptable  feedforward 
network.    It  will  be  found  that  a  number  of  acceptable  approximations 
can  be  formed  by  the  Pade'  method  to  meet  the  network  requirements. 
The  problem  will  be  to  choose  the  best  approximation  and  to  find  the 
conditions  for  which  it  is  the  best.    In  this  thesis,  a  frequency  domain 
technique  is  devised  for  choosing  the  best  approximation  in  a  minimum  inter 
gral    square  error  sense.    Experimental  examples  are  then  given  to  show 
the  effectiveness  of  the  approximations  and  the  compensation  method. 
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II.     FEEDFORWARD  COMPENSATION  OF 
SYSTEMS  WITH  TIME  DELAY 

Methods  of  Compensation.     Given  a  system  such  as  the  one  shown  in 
Fig.    1,  the  solution  to  the  stability  problem  involves  the  removal  of 
the  e-sl  term  from  the  denominator  of  the  closed  loop  transfer  function. 
Two  methods  of  doing  this  immediately  come  to  mind  and  are  shown  in 
Fig.  4.     The  transfer  function  for  system  1  (Fig.  4a)  can  be  expressed 
as 

C(s)   ^  G,(s)Gi.Me~slr 

It  can  be  seen  that  if  Gg(s)  is  set  equal  to  C2(l»e~sT)  the  transfer 
function  reduces  to 

■£&■  —  &(S)Ga(S)  _-sr     (4) 

It  is  recognized  that  this  equation  represents  nothing  more  than  the 
transfer  function  of  a  closed  loop  system  (without  a  delay  element  in 
the  loop)  cascaded  with  the  time  delay. 

In  system  2  (Fig.  4b),  if  G3(s)  is  again  set  equal  to  G2(l-e"sT), 
the  transfer  function  can  be  immediately  written  as 

C(s)    ___    ___Gjhlj^i)____ 

R(s)    -      /  v-  Gt(s)6a(s)  (5) 

It  can  be  seen  that  both  techniques  meet  the  objective  of  removing  the 
e   s     term  from  the  denominator  of  the  closed  loop  transfer  function. 


Both  systems  use  a  standard  network,  which  has  a  transfer  functi 


on 
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1-e        ,  cascaded  with  a  G?(s)  network.    This  means  that  the  time 
delay  network  needs  to  be  synthesized  only  once,  normalized  with 
respect  to  the  time  delay,  for  all  problems.     The  only  requirement  is  an 
impedance  match  between  the  time  delay  network,  G2(s)  network,  and 
any  other  sections  in  the  feedforward  path.     There  is  a  subtle  difference 
in  the  two  methods  of  compensation  however.    This  difference  involves 
the  nature  of  the  output.    The  two  systems  can  be  reduced  to  the 
equivalent  block  diagram  form  shown  in  Fig.  5.    The  output  of  the  first 
system  is  still  delayed  by  a  factor  T  while  the  output  of  the  second 
system  starts  to  respond  immediately  after  the  application  of  the  input 
signal.     First  impressions  would  indicate  that  the  second  system  is 
better  than  the  first  since  the  time  delay  has  been  completely  eliminated, 
However,  it  remains  to  be  seen  from  experimental  data  that  this  is  the 
case.    It  is  possible  that  noise  attenuation  characteristics,  speed  of 
response,  or  other  desirable  parameter   might  require' 
use  of  the  first  system.    A  practical  consideration  might  also  dictate 
the  use  of  the  first  system.    If,  for  example,  the  time  delayed  plant  is 
a  chemical  process,  it  would  be  impossible  to  "mix"  the  chemical 
output  with  anelectrical  feedforward  signal.     On  the  other  hand,  the 
electrical  feedforward  signal  could  be  added  to  the  electrical  feedback 
signal  coming  from  some  output  sensor,  as  required  in  the  first  system. 
The  Approximation  Requirement.    The  first  step  in  implementing  either 

of  these  two  systems  is  to  synthesize  a  network  which  generates  the 

—  sT 
transfer  function  1-e        ,  which  is  identical  to  the  transfer  function 

required  for  the  Smith  Linear  Predictor.     Since  it  is  a  transcendental 
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function,  some  form  of  approximation  must  be  used  if  the  transfer 
function  is  to  be  realized  by  lumped  passive  components.     Obviously- 
distributed  parameter  delay  lines  can  be  used  but  if  T  is  large,  as  can 
be  expected  in  a  process  control,  these  devices  are  not  practical.    A 
set  of  linear  filter  transfer  functions  can  be  derived  using  any  of  the 
standard  approximation  methods.    The  background,  procedure,  and 
accuracy  behind  the  Pade'  method  of  polynomial  approximation  is 
contained  in  the  next  chapter.    Approximations  based  on  this  method  are 
used  in  this  thesis  for  the  feedforward  network  required  in  the  two 
compensation  schemes  of  Fig.  4. 
Summary.    The  results  of  this  chapter  can  be  summarized  as  follows: 

1.  A  closed  loop  control  system  with  time  delay  can  be  stabilized 

if  the  proper  signal  can  be  fed  forward  around  the  time  delay  element, 

2.  The  filter  required  to  generate  the  feedforward  signal  has  the 
transfer  function 

G(s)  *  l-e"sT 

3.  Approximation  techniques  must  be  used  to  synthesize  the 
required  transfer  function  with  lumped,  linear  components. 


18 


III.     DESIGN  OF  OPTIMAL  FEEDFORWARD  FILTERS 

The  Pade'  Approximation.    A  required  step  in  the  synthesis  of  any 
function,  in  lumped,  linear,  passive,  time  invariant  form,  is  to  insure 
that  the  function  is  rational.    The  Pade'  approximation  is  a  method  of 

representing  a  polynomial  or  series  of  polynomials  as  a  rational 

3   4 
function.     '      Assume  that  some  function,  G(s) ,  is  to  be  approximated 

by  a  rational  fraction  of  an  n-th  order  polynomial  in  s,  divided  by  an 

m-th  order  polynomial  in  s.     G(s)  could,  for  example,  be  the  Taylor 

series  expansion  of  the  desired  transfer  function  l-e~s™ .    Thus  G(s) 

can  be  represented,  in  general  rational  form,  as 

',,        0.0  +  4,5+  '  •  -  +-a^sn 

G($)      —    ■■  ■  ■   "--~»-«»^ — — - — (6) 

fc0+  b,S+   •   •   •  +  hfoSm 

This  expression  may  be  expanded  by  cross  multiplication  to  give 

(both,  S  <h  •  ')(&0 1  Bt  &+>".) 

(7) 

where  Gn,G    ,  etc.  are  the  coefficients  of  the  Taylor  series  expansion 
of  G(s)  around  the  origin  in  the  s -plane.     By  equating  equal  powers  of 
s,  a  set  of  simultaneous  algebraic  equations  are  obtained.     For  a  given 
number  of  terms  in  the  Taylor  series  and  for  a  particular  number  of  terms 
in  the  rational  fraction,  there  is  a  finite  number  of  equations  which  are 
solved  simultaneously  to  obtain  the  a's  and  b's. 


19 


From  the  s     terms  Aq  =  bnGn 

s1  Al  =  blGo  +  bQGl 

s2  A2  =  b2GQ  +  biG],  +  b0G2 

s3  A3  =  b3G0  +  b2G2  +  bxG2  +  b0G3 


sn  An  =  bnG0  +  Vl^i  +  ' ' '  +  b0Gn  (8) 

The  rule  for  finding  the  sn  equation  should  be  obvious.    It  is  noted  that 
it  requires  m+n+1  terms  of  the  Taylor  series  to  find  n+1  (or  n)  coefficients 
for  the  numerator  and  m  (or  rrH-1)  coefficients  for  the  denominator  of  the 
rational  fraction,  and  that  m+n  simultaneous  equations  are  used  from  (8). 
This  still  allows  one  of  the  coefficients,  say  bg,  to  be  picked 
arbitrarily  and  the  only  difference  in  the  approximations  will  be  a 
constant  multiplier  for  different  values  of  b„.    It  is  realized  that  when 
the  rational  fraction  is  divided  out,  in  ascending  powers  of  s,  the  first 
m+n+1  terms  of  the  Taylor  series  are  reproduced  exactly.    Therefore, 
the  "maximum  error"  incurred  in  using  the  Pade'  approximation  is  no 
more  than  the  remainder  of  the  Taylor  series,  taken  after  m+n+1  terms. 
Nomenclature.    The  following  nomenclature  will  be  used  in  the  remainder 

of  this  thesis  to  refer  to  a  particular  approximation  or  transfer  function. 

-sT 

Emn  refers  to  the  Pade'  approximation,  of  the  function  e        ,  characterized 

by  an  n-th  order  polynomial  divided  by  an  m-th  order  polynomial  in  s. 

-sT 

But  e         is  not  the  transfer  function  required  in  the  feedforward  path  of 

the  systems  in  Fig.  4.    Therefore,  Gmn  will  refer  to  the  approximation 

— sT 

of  the  desired  transfer  function,   1-e        ,  derived  from  the  Pade1 
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approximation  of  e   D     of  the  same  order.     For  example,  if 


^21 

6 -'?::   '/■■" 

6<;~H$T*$sTa 

then 

(9) 


6  *  L/5T  +  S2T2 


(10) 


Derivation  of  the  Desired  Filter.    The  Taylor  series  expansion  of  e~x, 
about  the  value  x  =  0,  is  the  well  known  power  series 

e~*~  \-x+ztXz-ttXs-I""'  (11) 

This  is  the  series  that  is  used  to  derive  rational  approximations  for  the 

-  sT  ~x 

function  1-e        ,  with  sT  =  x.     Pade1  approximations,  for  the  function  e 

are  tabulated  in  many  sources.      Appendix  1(a)  gives  such  a  list  of 
approximations.    Appendix  1(b)  gives  the  derived  filter  transfer 
functions  derived  from  the  appropriate  approximation,  while  Appendix  1(c) 
gives  the  derived  filter  transfer  functions  in  factored  form,  which  is  not 
available  in  the  general  literature.    The  polynomials  were  factored  using 
the  POLYRT  subroutine  on  the  CDC  1604  digital  computer  and  rounded 
off  to  three  decimal  places.    While  the  factored  form  is  never  used 
directly  in  any  calculations  in  this  paper,  it  does  give  a  little  insight 
as  to  the  location  of  the  poles  and  zeros  of  each  of  the  filter  transfer 
functions.    The  poles  and  zeros  are  plotted  in  Fig.  6.    There  appears  to 
be  a  pattern  to  the  poles  and  zeros  as  the  order  of  s  in  either  the 
numerator  or  denominator  is  increased.    This  pattern  is  not  known  but 
if  it  were,  it  could  be  used  to  formulate  a  general  rule  for  finding  the 
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filter  transfer  function  of  any  desired  order.     However,  the  solution 
of  (8)  is  quite  straight  forward  in  solving  for  any  particular  Emn  so  that 
no  effort  was  expended  to  determine  an  expression  for  this  pole  zero 
pattern. 

There  is  no  guarantee  that  all  of  the  transfer  functions  found  in 
Appendix  1(b)  are  realizable  with  passive  elements.    The  properties  of 
a  transfer  function,  made  up  of  lumped,  linear,  passive  elements  are 
well  known  but  are  reiterated  here." 

1.  Poles  are  in  the  left  half  of  the  s  plane  or,  for  the  lossless 
case,  on  the  imaginary  axis  and  simple.    Complex  poles  must  be 
conjugate  pairs. 

2.  Zeros  may  have  any  s  plane  location  but  complex  zeros  must 
occur  in  conjugate  pairs.    A  minimum  phase  transfer  function  has 
all  of  its  zeros  in  the  left  half  of  the  s  plane. 

3.  Poles  may  not  be  located  at  the  origin  or  at  infinity. 

4.  The  degree  of  the  numerator  may  equal  but  not  exceed  the 
degree  of  the  denominator. 

Item  4  above  rules  out  all  transfer  functions  where  m  is  less  than  n. 
These  transfer  functions  will  no  longer  be  considered  in  this  thesis. 
Referring  back  to  Fig.   6,  it  is  seen  that  not  every  Pade'  derived  transfer 
function  can  be  synthesized  with  lumped,  linear,  passive  elements, 
even  though  m  is  greater  than  n.    As  the  order  of  s  in  the  denominator 
is  increased,  holding  the  order  of  s  in  the  numerator  constant,  the 
motion  of  the  poles  is  such  that  poles  begin  to  appear  in  the  right  half 
of  the  s  plane.    This  violates  rule  1  above  and  renders  these  transfer 
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functions  unacceptable  as  an  approximation  to  the  desired  filter 
transfer  function,   1-e        .    It  is  also  noted  that  for  every  transfer 
function  where  m  equals  n,  an  amplifier  is  required,  for  D.C.  gain,   in 
addition  to  the  passive  elements  for  synthesis. 

There  is  no  indication  as  to  which  of  the  filter  transfer  functions 
is  "the  best.  "    Intuitively,   G77  should  be  the  best  approximation  to  the 
ideal  transfer  function  since  its  derivation  required  more  terms  of  the 
power  series  than  any  of  the  others  in  the  list.     On  the  other  hand,  G?y 
requires  more  effort  to  synthesize  since  it  is  more  complex  than  the 
others.    It  may  be  that  a  simple  filter  will  work  just  as  well  .-as  a 
complex  one  in  the  feedforward  path  of  the  systems  in  Fig.  4.    A 
measurement  of  which  is  the  best  approximation  to  the  ideal  filter 
transfer  function,  1-e        ,  is  attempted  below. 

The  Minimum   Integral    Square   Error    Index.      The   index,    chosen   to  measure 
which  of   the   filter   transfer   functions      is    the  best   approximation,    is 

the  Minimum     Integral  Square  Value  of  the  difference  between  the 

-sT 

responses  of  the  ideal  transfer  function,   1-e        ,  and  the  Pade'  derived 

filter  approximation  G       .     The  integral  square  error  can  be  calculated   in 
either  the  time  or  frequency  domains.    Although  the  time  domain  has 
received  predominant  attention  in  the  literature,  several  advantages 
uncovered  during  this  investigation  urge  the  use  of  the  frequency  domain 
technique.     The  first,  and  most  important  of  the  advantages,  is  the 
ability  to  solve  for  a  frequency  spectrum  which  can  be  later  used  to  find   the 
integral    square  error  for  any  input.     This  spectrum,  called  the  weighting 
function,  gives  the  user  a  certain  feeling  for  the  "quality"  of  the 

25 


network  under  consideration.    When  time  domain  techniques  are  used, 
an  input  signal  must  immediately  be  chosen  and  the  result  applies  for 
that  input  only.    The  result  that  is  obtained  is  only  a  number  and  is 
not  a  very  satisfying  "answer. "    Without  other  networks  (numbers)  for 
comparison,  the  results  are  next  to  useless.    Even  with  several  networks, 
the  user  only  knows  which  one  has  the  smallest  integral  square  error,   not 
what  causes  that  error  or  what  effect  another  input  might  have.    This 
advantage  will  be  more  apparent  by  comparison  of  the  results  of  the 
two  techniques.    Another  advantage  of  the  frequency  domain  approach 
(or  disadvantage  of  the  time  domain  method)  is  not  having  to  take  the 
inverse  LaPlace  transformation  of  the  network  transfer  function.    Some 
of  the  larger  Pade'  approximations  are  very  difficult  and  time  consuming 
when  it  comes  to  obtaining  time  domain  solutions. 

In  view  of  the  above  arguments,  this  thesis  will  primarily  use  the 
frequency  domain  approach  to  selecting  the  best  Pade'  approximation 
for  use  as  a  filter  transfer  function  in  the  compensation  systems  of 

Fig.  4.    It  is  known  that  the  mean  square  value  of  a  signal  is  equal  to 

7 
the  value  of  its  autocorrelation  function  at  t  =  0. 

mean  sq.  value  of  v(t)  =        $\/y 


t^O 


§vAM  -lw 

(12) 

where  0vv(jw)  is  referred  to  as  the  Power  Density  Spectrum  of  the 
signal.    It  derives  this  name  from  the  fact  that  the  integral  in  (12)  equals 
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the  average  power  dissipated  in  a  1  ohm  resistor  by  the  signal 
voltage  v(t).   0vv(jw)  can  be  interpreted  as  a  lineardensity  distribution 
of  the  mean  square  value  of  the  signal  along  the  imaginary  axis  of  the 
s  plane,     y)  vv(jw)  is  obtained  from  the  relationship 

^(s)  s,T/fi      ^^  (13) 

and  then  setting  s  =  jw. 

Consider  the  hypothetical  experiment  shown  in  Fig,   7.     It  can  be 
shown    that  the  power  density  spectrum  of  the  error  signal  Y(s)  can  be 
expressed  as 


(&)  -      H(s)H(S)$„(s) 

-    His)  6 (~$)  $al(s) 

-i-  O(s)6(-S)0„d)  (i4) 


where 


Hfs)  "   l°e"ST  as) 


and  0]_j_(s)  is  l/2lltimes  the  LaPlace  transform  of  the  autocorrelation 
function  of  v^(t),  0  \2\s)  1S  l/2{jtimes  the  LaPlace  transform  of  the  cross 
correlation  function  of  v-^(t)  and  V2(t),  etc.    If  vj_(t)  equals  vpft),  as 
shown  by  the  dotted  lines  in  Fig.   7,  then  all  of  the  auto  and  cross 
correlation  functions  in  (14)  are  equal  and  will  be  designated  by  a 
single    0  vv(s)»     ^  *s  known  that  for  any  arbitrary  function  F 
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F (J"-1! '-i  ;■ 


F(f  UJ) 


Therefore,  from  (14)  with  s  =  jw 


(16) 


The  quantity  in  brackets  in  (17)  will  be  referred  to  as  W(jw) 


V 


(17) 


j>Wh    W(iw)   +  6(^-i!kJH(f^&^jj    (17a) 

By  trigonometric  substitution 

lam 

(\~CQS  UJT)   +f   SIM 


l_'3 


so  that 


(18) 


YHjuJ) 


(i- cas  ujt) 


(19) 


From  the  Pade'  approximations,  G       (s)  is  of  the  form 


ftn5<tAn-iS    "V  *"+/)!&  + Ao 


/;?<-* 


where  m  and  n  depend  on  the  particular  approximation  being 


considered.     Or  in  simpler  terms 


(B  (s) 


M*  (s)  +  /Vfi  ($) 
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(20) 


(21) 


where  the  M's  refer  to  the  sum  of  the  terms  with  even  powers  of  s  and 
the  N's  refer  to  the  sum  of  the  terms  with  odd  powers  of  s.     For  real 
frequencies  (when  s  =  jw) 

Gfyw)=      —L——J-—  =  fHtunfbfa)        (22) 


and 


a~         fiB    *    Ba  (23) 


(24) 


<S(rfWj 

From  (18)  and  (22),  it  follows  that 

-fiMJh  Cos  wf}  +Blw)[s>tii  wfj 

Substituting  (19),  (23)  and  (24)  into  the  bracketed  portion  of  (17)  yields 

"Z jfilW)  (l-cos  k»t)  +B(»J)Stff  t'jfj 

In  shortened  form,   (17)  becomes 

§yr(AfJ)  =  Wfc&)$n,tfwJ  (26) 

W(jw)  is  given  the  name    Spectral  Weighting  Function.     Referring  to 
Fig.  7,  the  spectral  weighting  function  does  just  what  its  name  implies. 
It  weights  the  magnitude  of  the  power  density  spectrum  of  some  input 
signal  v(t)  to  yield  the  power  density  spectrum  of  the  error.    All  that  is 
now  required,  to  determine  which  of  the  Pade'  derived  filter  transfer 
functions  is  the  best,  is  the  substitution  of  the  various  approximations, 
Gmn,  into  (25)  to  determine  W(jw) .     Multiplying  the  spectral  weighting 
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function  by  the  power  density  spectrum  of  an  appropriate  input  signal 
yields  the  power  density  spectrum  of  the  error  for  that  particular 
transfer  and  input  function.     Integrating  the  power  density  spectrum 
according  to  (12)  yields  the  integral   square  value  of   the  error   function. 
The  Pade'  derived  filter  transfer  function  that  has  the  smallest  area 
under  the  spectrum,   i.e.   the  minimum  integral   square  error,    is   the  one 
which  is  the  best  approximation  to  the  ideal  transfer  function,   1-e        . 
The  Spectral  Weighting  Function.     In  order  to  generalize  the  selection 
process  for  any  input  signal,  the  spectral  weighting  functions,  for  each 
of  the  acceptable  Pade'  derived  transfer  functions,  should  be  presented 
in  graphical  form.     This  permits  the  user  to  know  in  advance  which 
frequencies  of  the  input  signal  the  various  transfer  functions  weight 
the  most.    The  transfer  function,  that  best  suppresses  the  frequencies 
in  the  input  signal,   is  the  one  that  will  have  the  minimum  integral  square 
error.     It  is  significant  to  note  that  the  best  approximation  depends 
upon  the  spectral  distribution  of  the  input  signal  as  well  as  the  spectral 
weighting  function.   Due  to  the  large  number  of  Pade'  derived  filter 
transfer  functions  and  the  complexity  of  the  spectral  weighting  function 
equation,  these  curves  were  computed  using  a  digital  computer.    A 
program  was  written  to  compute  and  plot  W(jw)  and  is  presented  in 
Appendix  2(b).    The  results  are  shown  in  Appendix  3.    The  horizontal 
axis  for  the  graphs  in  Appendix  3  is  plotted  as  the  log-^ q  of  a  normalized 
frequency  time  delay  product,  i.e.  x  =  log  to  wT.     In  general,  the 
spectral  weighting   functions   for   the  Pade1    derived   filter   transfer 
functions  have  a  larger  weighting   effect   at  higher   frequencies   and 
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appear    to  act  like  high  pass  filters  of  cutoff  frequency,  w«.    With 
the  exception  of  the  curves  corresponding  to  the  case  when  m  equals  n, 
all  of  the  spectral  weighting  functions  settle  to  a  constant  amplitude  of 
one.    The  spectral  weighting  function  of  a  transfer  function  having  m 
equal  to  n  has  a  constant  amplitude  oscillation  about  the  value  two. 
If  the  frequency  components  of  an  input  signal  were  lower  than  some  Wq, 
where  the  spectral  weighing  function  is  zero,  then  the  integral  square 
error  would  be  zero.     For  such  an  input  signal,  the  Pade'  derived  filter 
transfer  functions  would  be  an  excellent  approximation  to  the  ideal 
transfer  function. 

Comparison  of  Filter  Transfer  Functions.    The  input  signal  that  is  quite 
often  used  to  test  the  performance  of  a  control  system  or  filter  is  the 
unit  step.    The  idea  is  that  any  signal  can  be  approximated  by  a  sum 
of  little  steps.    Any  filter  or  system  that  would  perform  well  with  a 
single  step  input  should  work  for  a  sum  of  steps.    A  unit  step  may  not, 
however,  be  a  very  good  choice  for  an  input  signal  in  the  selection 
process  described  above.     Note  that  the  minimum  integral  square  error  is 
a  nonlinear  performance  index  and  the  principle  of  superposition  cannot 
be  used.    The  Pade1  derived  filter  transfer  function  that  is  selected  as 
the  best,  based  on  a  unit  step,  cannot  be  assumed  to  be  the  best  for 
some  Other  input.     The  minimum  integral  square  error  performance   index 
does  not,  however,  rule  out  the  possibility  that  the  best  transfer 
function  for  a  step  input  may  also  be  best  for  other  inputs.    It  all  depends 
on  the  power  density  spectrum  of  the  input.     For  the  purpose  of  this  paper, 
a  unit  step  will  be  used  as  the  input  signal  when  selecting  a  Pade' 
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approximation  to  use  as  the  feedforward  filter  transfer  function.    This 
choice  was  made  simply  because  the  unit  step  is  general  and  easy  to 
work  with.    In  an  actual  problem,  where  the  input  signal  is  known,  the 
power  density  spectrum  of  the  actual  input  would  be  used.    Also,  the 
unit  step  makes  the  examination  of  experimental  results  a  familiar  task, 
The  power  density  spectrum  of  a  unit  step  is  given  by  the  equation 

(i)        \  •  '  '    '  <■'      '  (27) 


The  digital  computer  is  also  used  to  multiply  the  spectral  weighting 
function  by  the  power  density  spectrum  of  a  unit  step  and  to  numerically 
integrate  the  resulting  area  under  the  power  density  spectrum  of  the 
error.    The  program  that  was  written  to  perform  this  operation  is  given  in 
appendix  2(c).    The  integral    square  error  for  each  of  the  acceptable  Pade' 
derived  transfer  functions  is  tabulated  in  table  1.    The  table  is  broken 
into  several  bands  to  show  how  the  frequencies  of  the  spectrum  cont- 
ribute to  the  integral  square  error.     Figs.  8  and  9  are  graphical 
representations  of  the  data  in  table  1.    The  solid  lines  in  Fig.   8  show  the 
variation  in  the  integral  square  error  as  one  of  the  indices  (m  or  n)  is 
varied.    Horizontal  movement  along  one  of  the  solid  lines  from  one 
point  to  another  represents  a  change  in  the  order  of  s  in  the  numerator, 
while  holding  the  order  of  s  in  the  denominator  constant.     It  is  observed 
that,  as  long  as  n  is  less  than  rn,  an  increase  in  the  order  of  s  in  the 
numerator  results  in  a  decrease  in  the  integral  square  error.    The  amount 
of  the  decrease  becomes  smaller  as  n  approaches  m.    When  n  becomes 
equal  to  rn ,  there  is  a  substantial  increase  in  the     integral  square  error  and 
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Table  1 

Integral  Square  Errors  (Frequency  Domain)  for 

Pade'  Derived  Filter  Transfer  Functions 
Step  Input  -  Time  Delay    1  Second 


Band  0        less  than  .01  rps  Band  3  10.  to  100.  rps 

Band  1         .01  to  1.  rps  Band  4  100.  to  1000.  rps 

Band  2         1.  to  10.  rps  Band  5  over  1000.  rps 


Band  0  12  3  4  5  M.S.E. 


G  10  aprx.O  .0165 

G  11  aprx.O  .0003 

G  20  aprx.O  .0016 

G  21  aprx.O  .0000 

G  22  aprx.  0  aprx.O 

G  30  aprx.O  aprx.O 

G  31  aprx.O  aprx.O 

G  32  aprx.  0  aprx.O 

G  33  aprx.O  aprx.O 

G  40  aprx.O  aprx.O 

G  41  aprx.  0  aprx.  0 

G  42  aprx.  0  aprx.  0 

G  43  aprx.O  aprx.O 

G  44  aprx.O  aprx.O 

G  51  aprx.  0  aprx.  0 

G  52  aprx.O  aprx.O 

G  53  aprx.  0  aprx.O 

G  54  aprx.  0  aprx.  0 

G  55  aprx,  0  aprx.  0 

G  61  aprx.  0  aprx.  0 

G  62  aprx.  0  aprx.O 

G  63  aprx.  0  aprx.O 

G  64  aprx.O  aprx.O 

G  65  aprx,  0  aprx.O 

G  66  aprx.O  aprx.O 

G  72  aprx.O  aprx.O 

G  73  aprx.O  aprx.O 

G  74  aprx.O  aprx.O 

G  75  aprx.O  aprx.  0 

G  76  aprx.O  aprx.O 

G  77  aprx.O  aprx.O 


.1868 

.0293 

.0029 

aprx.  0 

.2355 

.2021 

.0717 

.0057 

aprx.  0 

.2699 

.1142 

.0286 

.0029 

aprx .  0 

.1471 

.0740 

.0203 

.0029 

aprx.  0 

.1060 

.0974 

.0506 

.0057 

aprx.  0 

.1537 

.0974 

.0286 

.0029 

aprx,  0 

.1289 

.0449 

.0283 

.0029 

aprx.  0 

.0761 

.0369 

.0290 

.0029 

aprx,  0 

.0687 

.0401 

.0605 

.0057 

aprx.  0 

.1063 

.1474 

.0286 

.0029 

aprx.  0 

.1789 

.0347 

.0286 

.0029 

aprx.O 

.0662 

.0197 

.0296 

.0029 

aprx,  0 

.0522 

.0137 

.0340 

.0029 

aprx.  0 

.0500 

.0095 

.0695 

.0057 

aprx.O 

.0009 

.0366 

.0288 

.0029 

aprx.  0 

.0682 

.0138 

.0292 

.0029 

aprx.O 

.0450 

.0064 

.0306 

.0029 

aprx.  0 

.0398 

.0030 

.0344 

.0029 

aprx.  0 

.0402 

.0012 

.0584 

.0057 

aprx.  0 

.0652 

.1366 

.0286 

.0029 

aprx.O 

.1681 

.0130 

.0286 

.0029 

aprx.O 

.0446 

.0040 

.0204 

.0029 

aprx.  0 

.0353 

.0012 

.0292 

.0029 

aprx.  0 

.0322 

.0002 

.0300 

.0029 

aprx,  0 

.0332 

.0000 

.0486 

.0057 

aprx,  0 

.0544 

.0234 

.0284 

.0029 

aprx.O 

.0545 

.0036 

.0272 

.0029 

aprx.O 

.0338 

.0006 

.0252 

.0029 

aprx.O 

.0207 

.0000 

.0242 

.0029 

aprx.O 

.0271 

aprx.  0 

.0254 

.0029 

aprx.O 

.0203 

aprx.Q 

.0410 

.0057 

aprx.O 

.0467 
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when  n  becomes  larger  than  m  the  integral   square   error  goes   to   infinity. 
Moving  vertically  from  one  solid  line  to  another  represents  a  change 
in  the  order  of  s  in  the  denominator  while  holding  the  order  of  s  in  the 
numerator  constant.     For  these  changes,  the  integral  square  error  usually 
decreases  as  the  order  of  s  is  increased.    The  exceptions  to  this 
occur    when  there  is  a  pole  "close"  to  the  imaginary  axis  (Fig.   6)  such 
as  G4Q  and  Ggj .    The  dotted  lines  in  Fig.  8  represent    the  change  in  the 
integral    square  error  as  both  indicies  are  changed,  but  the  difference 
is  held  constant.    It  is  readily  observed  that  the  amount  of  improvement 
in   the   integral   square  error   is   exponentially  decreasing  as   the  order    (and 
complexity)  of  the  transfer  function  is  increasing.    Figs.  9  (a  and  b) 
show  the  changes  in  the  contribution  to  the  integral  square  error  by  the 
various  frequency  bands.    It  is  observed,  from  Fig.  9(a),  that  changing 
the  order  of  s  in  the  denominator,  while  holding  the  order  of  s  in  the 
numerator  constant,  has  little  effect  on  which  band  contributes  the  most 
to  the  integral     square  error.     On  the  other  hand,  Fig.   9(b)  shows 
that  changing  n  while  holding  m  constant  has  a  considerable  effect  on 
the  frequency  band  contributing  the  most  area.    As  n  is  increased,  the 
higher  frequency  bands  tend  to  contribute  more.    This  fact  could  be  used 
to  advantage  in  many  systems.    A  large  number  of  closed  loop  systems 
tend  to  act  as  a  low  pass  filter.    If  one  of  the  Pade'  derived  filters  is 
used  in  such  a  closed  loop  system,  the  low  pass  filtering  effect  of  the 
system  should  suppress  the  higher  frequencies  that  are  weighted  more 
by  the  high  pass  feedforward  filter.    This  means  that  the  integral  square 
error  for  a  unit  step  input  to  such  a  system  would  be  less  than  that 
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when  n  becomes  larger  than  m  the  integral   square   error   goes   to   infinity. 
Moving  vertically  from  one  solid  line  to  another  represents  a  change 
in  the  order  of  s  in  the  denominator  while  holding  the  order  of  s  in  the 
numerator  constant.     For  these  changes,  the  integral  square  error  usually 
decreases  as  the  order  of  s  is  increased.    The  exceptions  to  this 
occur    when  there  is  a  pole  "close"  to  the  imaginary  axis  (Fig.   6)  such 
as  G4Q  and  Gg^.    The  dotted  lines  in  Fig.  8  represent    the  change  in  the 
integral    square  error  as  both  indicies  are  changed,  but  the  difference 
is  held  constant.    It  is  readily  observed  that  the  amount  of  improvement 
in   the   integral  square  error   is   exponentially  decreasing  as    the  order    (and 
complexity)  of  the  transfer  function  is  increasing.     Figs.  9  (a  and  b) 
show  the  changes  in  the  contribution  to  the  integral  square  error  by  the 
various  frequency  bands.    It  is  observed,  from  Fig.   9(a),  that  changing 
the  order  of  s  in  the  denominator,  while  holding  the  order  of  s  in  the 
numerator  constant,  has  little  effect  on  which  band  contributes  the  most 
to  the  integral     square  error.     On  the  other  hand,  Fig.  9(b)  shows 
that  changing  n  while  holding  m  constant  has  a  considerable  effect  on 
the  frequency  band  contributing  the  most  area.    As  n  is  increased,  the 
higher  frequency  bands  tend  to  contribute  more.    This  fact  could  be  used 
to  advantage  in  many  systems.    A  large  number  of  closed  loop  systems 
tend  to  act  as  a  low  pass  filter.    If  one  of  the  Pade'  derived  filters  is 
used  in  such  a  closed  loop  system,  the  low  pass  filtering  effect  of  the 
system  should  suppress  the  higher  frequencies  that  are  weighted  more 
by  the  high  pass  feedforward  filter.    This  means  that  the  integral  square 
error  for  a  unit  step  input  to  such  a  system  would  be  less  than  that 
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calculated  above.  It  would  therefore  be  advantageous  to  use  the  filter 
transfer  function  with  the  highest  frequency  weighting  in  systems  that 
act  as  low  pass  filters. 

The  Integral  Square  Error  Index  .(Time  Domain).    In  order   to  cheek  the 
above  results,    the   integral   square   error   is   calculated   for   a   few  cases 
in   the   time   domain.      The   same  hypothetical   experiment,    shown   in  Fig.    7, 
is  used.      The   equation  for   the   error"  signal   is 
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Taking  the  inverse  La  Place  Transform 
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The  integral  square  error  is  given  by  the  equation 
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(32) 


It  is  easily  seen  that  the  first  integral  reduces  to  T  (the  area  under  a 
pulse  of  magnitude  1  and  duration  T) .    The  second  integral  is  equal  to 
the  area  under  the  response  curve  of  the  Pade'  derived  transfer  function, 
Gmn,  from  time  0  to  T.    The  last  integral  represents  the  area  under  the 
square  of  the  filter  response  curve  for  all  time.    The  lower  limits  can 
be  changed  to  0  since  the  response  does  not  start  prior  to  the  application 
of  the  input  at  time  0.    The  equation  for  the  integral  square  error,  for  a 
unit  step  input,  therefore  reduces  to 


J    [         I     e'  c-  I  •   \f\  TV'      /     /,  *G/H    V,  .  (33) 
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The  digital  computer  program  presented  in  appendix  2(d)  was 

written  to  compute  I  from  (33).     Because  of  the  long  computation  time, 
only  a  few  of  the  Pade'  derived  filter  transfer  functions  were  tested  by 
this  method.    The  results  are  given  in  table  2.    Since  all  of  the 
approximations  contain  decaying  exponential  terms,  the  upper  limit  for 
the  third  integral  in  (33)  was  set  at  20  times  the  longest  time  constant. 
By  comparing  tables  1  and  2,  it  is  seen  that  the  results  of  the  two 
methods  agree  to  within  the  third  decimal  place.     Better  agreement  can 
be  obtained  if  the  tolerance  on  the  numerical  integration  is  made  smaller. 

Such  action  would  increase  the  computation  time  however. 

-sT 

Other  Approximations  for  1-e  Several  other  sets  of  approximations 

were  investigated  but  none  performed  as  well  as  the  Pade'.    The  defining 

o 

equations  for  the  other  approximations  investigated  are 

/ 
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Table  2 


G  10 
G  11 
G  20 
G  21 
G  22 


Integral  .Square  Errors    (Time  Domain)    for   Several 
Pade'  Transfer  Functions  -  Step  Input 
Time  Delay  -  1  Second 


first 
integral 


second 

third 

inteqral 

integral 

.6321 

.5000 

.8647 

1.0000 

.8006 

.7490 

.8818 

.8696 

.9225 

.9992 

J.S.E, 

.2358 
.2707 
.1477 
.1059 
.1543 


Table  3 


Integral  Square  Errors  for 


Bj  to  B     and  C     to  C5 
Step  Input  -  Time  Delay  1  Second 


Band  designations  are  the  same  as  for  Table  1 


band 


M.S.E. 


B  1 

aprx.  0 

.0166 

.1868 

.0293 

.0029 

aprx.  0 

.2355 

B  2 

aprx.  0 

.0054 

.1289 

.028  6 

.0029 

aprx.  0 

.1662 

B  3 

aprx.  0 

.0025 

.1012 

.0284 

.0029 

aprx.  0 

.1353 

B  4 

aprx.  0 

.0016 

.0840 

.0286 

.0029 

aprx.  0 

.1168 

B  5 

aprx.  0 

.0010 

.0719 

.0286 

.0029 

aprx.  0 

.1044 

C  1 

aprx.  0 

,1248 

.1808 

.0608 

.0057 

aprx.  0 

.4676 

C  2 

aprx.  0 

.2779 

.2852 

.0512 

.0057 

aprx.  0 

.6201 

C  3 

aprx.  0 

.2928 

.42  62 

.0633 

.0057 

aprx.  0 

.  7884 

C  4 

aprx.  0 

.2986 

.5621 

.0589 

.005  7 

aprx .  0 

.9253 

C  5 

aprx.  0 

.3014 

.6621 

.0484 

.0057 

aprx.  0 

1.0176 
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A 
These  were  chosen  mainly  on  the  basis  that  they  are  relatively  easy  to 

synthesize.     For  example,   (35)  is  a  lattice  network  yielding  the  proper 
time  delay.    The  results  of  the  investigation  are  given,  in  table  3.    The 
format  is  the  same  as  that  for  table  1  (i.e.  the  frequencies  in  each 
band  are  the  same  as  for  the  bands  in  table  1) .    It  is  noted  that  the  integ- 
ral square    error  for  the  B    approximation  becomes  smaller  as  n  is 
increased.    However,  for  the  same  degree  of  complexity,  the  Pade* 
approximation  has  a  smaller  mean  square  error.    The  Cn  approximation 
does  not  perform  well  at  all.     The  smallest  integral  square  error  occurs 
when  n  =  1  and  as  n  is  increased,  the  integral  square  error  increases.      The 


smallest   integral   square   error   for  a   C     approximation   is   larger   than   the 
largest   integral   square   error   for   a  Pade'    approximation. 
Summary.    From  the  above  observations,  it  appears  that  the  Pade' 
approximation  derived  from  the  largest  number  of  terms  in  the  power 

series  has  the  smallest   integral   square   error   and   is   therefore   the  best 

-sT 
approximation  to  the  ideal  transfer  function,   1-e        .    This  statement 

is,  of  course,  subject -to  the  condition  that  the  denominator  is  of  higher 
order  than  the  numerator  and  that  there  are  no  poles  in  the  right  half  of 
the  s  plane.    It  is  speculated  however,  that  the  use  of  this  approximation 
may  not  be  the  most  practical  in  many  cases.    The  increase  in  complexity 
(and  therefore  in  difficulty  in  synthesis)  may  far  outweigh  the  advantage 
Of  a  small  decrease  in  the    integral    square   error.      As   stated   before,    the 
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selection  of  how  complex  a  Fade1  derived  filter  transfer  function  to  use 
will  depend  on  the  system  being  compensated,  the  duration  of  the  time 
delay,  the  nature  of  the  input  signal,  etc.    The  results  of  this  chapter 
can  be  summarized  as  follows: 

1.  The  Pade1  approximation  is  a  method  for  representing  a 
transcendental  function  as  a  rational  fraction  of  polynomials.    It 
is  found  by  the  simultaneous  solution  of  a  series  of  algebraic 

equations. 

-sT 

2.  Some  of  the  Pade1  approximations  of  the  desired  function  1-e 

can  be  realized  with  lumped,  linear  passive  elements. 

3.  The  best  approach  in  choosing  a  filter  transfer  function  is' 
through  the  frequency  domain.    In  this  procedure,  a  quantity  called 


the  Spectral  Weighting  Function   is   generated   from   the  desired   transfer 
function  and  the  transfer  function  of  a  filter  under  consideration. 
The  spectral  weighting   function  graphically  shows  which   frequencies 
in  the  input  signal  contribute  most  to  the   integral  square  error  of   the 
filter.    When  the  spectral  weighting    function  is  multiplied  by  the 
power  density  spectrum  of  an  input  signal,  the  power  density 
spectrum  of  the  error  is  obtained.    The  filter  with  the  minimum  area 
under  the  power  density  spectrum  of  the  error  is  the  one  with  the 
minimum   integral   square   error   and   is   considered   "the  best." 
4'.     The  Pade'  approximations  have  the  smallest  integral  square  error 
of  all  of  the  approximations  investigated  for  the  step  input. 
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IV.     EXPERIMENTAL  INVESTIGATIONS 

In  order  to  examine  the  effects  of  approximation  in  the  feedforward 
compensation  schemes  outlined  earlier,  a  hybrid  computer  model  was 
tested  with  several  of  the  Pade1  derived  filters  in  the  feedforward  path. 
The  specific  objectives  of  the  tests  were  to  ascertain  the  derogatory 

effects  of  using  an  approximation,  rather  than  the  ideal  transfer  function 

~sT 
1-e        ,  and  to  determine  what  effects  the  order  of  the  approximation 

has  on  the  time  domain  performance  of  the  feedforward  filter. 

The  Plant.     For  simplicity,  the  plant  described  in  the  introduction  was 

the  one  used  for  the  hybrid  computer  model.    The  plant  is  shown  again 

in  Fig.   10(a).     Fig.   10(b)  shows  the  computer  model.    An  EAC  model 

TR-20  analog  computer  was  used  to  simulate  the  linear  portions  of  the 

plant.    A  CDC  model  160  digital  computer  was  used  as  the  plant  time 

delay  element.    In  addition,  a  portion  of  the  analog  computer  provided 

signal  conditioning  for  proper  A/D  and  D/A  conversion.    The  A/D 

converter  requires  all  signal  inputs  to  be  between  0  and  -10  volts.    The 

analog  computer  setup  is  shown  in  Fig.   11.    The  maximum  allowable 

signal  at  any  point  in  the  computer  model  is  chosen  to  be  +8  volts. 

This  is  a  convenient  choice  of  limits  since  the  digital  computer 

operates  on  the  octal  number  system  and  a  direct  correlation  between 

Voltage  levels  in  the  analog  computer  and  numbers  in  the  digital 

computer  memory  results. 

The  Digital  Time  Delay.    The  basic  idea  behind  the  digital  time  delay 

is  not  complex.    The  analog  computer  output  is  sampled  periodically 
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and  the  number  corresponding  to  the  amplitude  of  that  sample  is  stored 
in  the  computer  memory.    At  the  end  of  the  time  delay  period  T,  the 
number  is  returned  to  the  analog  computer  as  a  voltage  signal  through 
a  hold  element.    If  the  sampling  rate  is  fast  enough,  the  voltage 
returned  to  the  analog  computer  appears  to  be  a  continuous  time  varying 
signal.    A  flow  graph  of  this  operation  is  shown  in  Fig.   12.     The 
computer  program  written  to  perform  this  operation  is  given  in  appendix  2(f) 
There  is  a  difficulty  in  the  use  of  the  program  listed  in  appendix  2(f) 
however.    As  the  time  delay  is  increased,  the  sampling  rate  is  decreased. 
This  program  is  only  adequate  for  for  delays  of  several  minutes  or  less. 
For  time  delays  longer  than  this,  the  output  of  the  D/A  converter  is  better 
described  as  piecewise  constant  rather  than  continuous. 
Simulation  of  a  Typical  Pade'  Derived  Filter.    The  most  difficult  step  in 
setting  up  a  computer  model  of  the  system  is  the  simulation  of  the 
feedforward  filter,  G      .     Large  gains  (larger  than  10'  in  G77)  are 
required  in  all  but  the  simplest  of  the  Pade'  derived  filters.    The  steps 
taken  to  simulate  G4~  are  given  here  as  an  example  of  the  procedure. 

_         _    8ff05T^aO{ST}2"fa,0(ST}3^(sT^ 

V3       8L!0-tH8OST+l£O(ST)s-hlG>(sT)\($T)'l 

Dividing  the  numerator  and  denominator  of  (36)  by  (sT)     and  setting 
equal  to  the  ratio  of  the  output  to  the  input. 


Go         8iO/^&6O7^*£0?sT$  +  I  (37) 

™%7?     'u°WflW(STf '  lb5T     l 
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By  cross  multiplication  and  rearranging  terms 


LA 


(38) 


<     ■"     -.-■■'   / 


(38)  is  in  the  proper  form  for  the  simulation  on  the  analog  computer. 
The  analog  computer  setup  for  G43  is  shown  in  Fig.   13.    The  other  Pade' 
derived  filter  approximations  are  simulated  in  the  same  way. 
The  Step  Response.    It  was  stated  in  the  introduction  that  the  system 
under  consideration  was  stable  with  a  time  delay  of  1  second  but 
unstable  if  the  time  delay  was  increased  to  2  seconds.    The  Bode 
diagram  for  the  plant  was  presented  in  Fig.   2  to  support  this.    To  verify 
this,  the  feedforward  filter  was  disconnected  in  the  computer  model  and 
the  model  was  subjected  to  a  1  volt  step.    The  results  are  shown  in 
Fig.   14.    The  plant  is  definitely  unstable  for  a  time  delay  of  2  seconds. 

To  meet  the  objectives  of  the  tests,  three  of  the  Pade1  derived 
filter  transfer  functions  were  used  as  the  feedforward  network  in  the 
computer  model.    The  three  that  were  used  are  G21/  G43,  and  G55. 

These  were  chosen  because  they  had  the  minimum  integral  square  error 
for  a  given  value  of  m.     The  time  delay  was  set  to  2,   3,  and  4  seconds. 
At  all  three  of  the  time  delay  settings,  the  plant  without  compensation 
would  be  unstable.    The  results  of  the  tests,  for  both  compensation 
methods,  are  shown  in  Figs.    15  (a,  b  and  c) .     It  is  observed  that  all 
three  Pade'  derived  filter  transfer  functions  stabilize  the  system  for  all 
three  values  of  time  delay  tested.    For  the  longer  delays,  the  more 
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complex  filters  do  a  better  job  with  respect  to  the  shape  of  the  response 
curve.    The  response  with  the  Gfit-  filter  and  with  a  4  second  time  delay- 
looks  just  like  the  response  of  the  uncompensated  system  with  no  time 
delay,  with  the  exception  of  the  time  shift.    It  is  interesting  to  note  that 
the  two  compensation  methods  (Fig.  4)  do  not  give  the  same  compensated 
system  output.    The  best  example  of  observing  this  is  the  response  of 
G?,  with  a  4  second  delay,  Fig.   15(a).    The  response  due  to  compensa- 
tion method  2  oscillates  much  more  than  the  response  due  to  compensa- 
tion method  1  and  the  settling  time  is  longer.    On  the  whole,  it  appears 
that  compensation  method  1  performs  better  than  compensation  method  2. 
This  is  very  fortunate  since  it  was  noted  earlier  that  practical 
considerations  might  require  the  use  of  this  compensation  scheme. 
Sensitivity.     No  extensive  investigation  was  made  into  the  area  of 
sensitivity.     Only  one  short  experiment  was  performed  to  get  an  idea  of 
the  tolerance  limits  on  the  filter  time  delay  setting.    The  plant  was 
simulated  on  the  hybrid  computer  with  a  two  second  time  delay.    This 
is  an  unstable  condition.    The  plant  was  then  compensated,  using 
compensation  method  1  and  the  Fade'  derived  filter  Got.    The  filter  time 
delay  was  varied  from  1  to  5  seconds,  in  half  second  steps.    The  1 
second  setting  was  the  shortest  delay  that  could  be  conveniently  set 
on  the  computer  and  at  5  seconds,  the  system  was  at  the  margin  of 
stability.    The  results  of  some  of  the  other  trials  are  shown  in  Fig.   16. 
It  is  observed  that  the  system  is  stable  for  a  large  range  of  settings  of 
filter  time  delay.    The  quality  of  the  response  deteriorates  as  the 
difference  between  the  plant  and  filter  time  delays  increases  but  the 
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V.     CONCLUSION 

It  was  easily  shown,  in  Chapter  2  of  this  thesis,  that  the 
feedforward  compensation  techniques  presented  in  Fig.  4  would  remove 
the  phase  lag  due  to  a  time  delay  in  a  closed  loop  control  system  and 
would  therefore  improve  the  stability  of  such  a  system.    It  was  then 
implied  that  the  required  feedforward  network  could  not  be  synthesized 
with  lumped,  linear,  passive,  time  invarient  circuit  elements.    The 
problem  was  therefore  reduced  to  one  of  finding  a  rational  approximation 
that  would  perform  as  near  to  the  unrealizable  ideal  circuit  as  possible. 
Since  the  function  to  be  synthesized  is  a  common  one,  a  number  of 
approximations  and  approximation  techniques  were  available  to  choose 
from.    The  performance  index,  chosen  to  rate  the  various  approximations, 
was  the  minimum  integral   square   error  between   the   step   response   of   the 
ideal  and  approximation  filters.    The  Pade'  method  of  polynomial 
approximation  was  finally  used  since  it  was  simple  to  apply  and  it  gave 
a  wide  variation  in  transfer  function  complexity.    This  large  variety  of 
circuits  means  that  a  large  range  of   integral  square  error  was  available  to 
test  the  effects  of  approximation  in  the  two  compensation  schemes.    It 
was  found,  through  computer  simulations,  that  approximations  of  the 
ideal  transfer  function  effectively  stabilize  systems  with  time  delay. 
Even  poor  approximations,  i.e.  filters  with  large  integral  square  errors, 
were  able  to  effect  a  degree  of  stability  in  such  systems.    In  addition, 
it  was  found  that  the  stabilizing  properties  of  the  two  compensation 
schemes  were  not  affected  by  small  changes  in  the  time  delay  setting 
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system  remains  stable.    This  feature  is  useful  in  situations  where  the 
plant  time  delay  is  either  not  known  exactly  or  it  changes. 
Summary.    The  observations  made  in  this  chapter  can  be  summarized 
as  follows: 

1.  The  Pade'  derived  filter  transfer  functions  make  effective 
feedforward  filters  in  the  compensation  schemes  outlined  in  Fig.  4. 
As  anticipated,  the  more  complex  transfer  functions,  which  are 
better  approximations  to  the  ideal  transfer  function,  give  a  better 
system  response. 

2.  Compensation  method  1  (Fig.  4(a))  gives  a  better  system  response 
than  compensation  method  2. 

3.  The  ability  of  the  Pade1  derived  filter  transfer  functions  to 
stabilize  an  unstable  plant  is  not  very  sensitive  to  changes  in  the 
filter  time  delay  setting. 
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of  the  approximation  feedforward  filter.    In  general,  it  can  be  concluded 
that  both  the  compensation  methods  presented  in  Chapter  2  are  excellent 
compensators  for  time  delay  and  that  the  required  feedforward  circuits 
are  not  very  sensitive  to  approximation  effects. 

The  concept  of  a  frequency  domain  analysis  of  the  minimum  integral 
square  error  criterion,  as  developed  in  Chapter  3,  is  as  important  as 
the  idea  of  feedforward  compensation.    In  most  of  the  literature,  the 
time  domain  analysis  is  used.    While  the  time  domain  technique  gives 
satisfactory  answers,  i.e.   gives  the  integral  square  error  for  a  problem 
under  consideration,  the  frequency  domain  analysis  gives  this 
information  and  much  more.    The  spectral  weighting  function,   an  inter- 
mediate step  in  the  frequency  domain  analysis,  gives  the  user  a 
graphical  indication  of  the  frequencies  that  are  contributing  to  integral 
square  error.    This  means  that  the  user  can  find  a  solution  for  his 
problem  based  on  the  frequency  spectrum  of  the  input  signal  and  the 
spectral  weighting  functions  of  the  various  solutions  under  consideration. 
The  spectral  weighting  function  can  be  computed  once  and  used  to 
calculate  the    integral   square   error   for   any   input.      These  advantages   are 

not  available  when  the  time  domain  analysis  is  used. 

Topics  for  Further  Study.    The  ideas  used  in  this  thesis  can  be  projected 
into  additional  areas  of  study.    For  example,  work  can  be  done  in  the 
case  where  the  time  delay  cannot  be  isolated  from  the  remainder  of  the 
open  loop  transfer  function  of  the  plant.    Also,  the  idea  of  a  sensitivity 
investigation  could  receive  additional  attention.    It  is  interesting  to  note 
that  the  ideal  feedforward  transfer  function  is  very  similar  to  the 
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transfer  function  of  a  Zero  Order  Hold  circuit,   (l™e"s-^)/s.    A  number 
of  the  Fade'  derived  filters  found  in  appendix  1  could  possibly  be  used 
in  conjunction  with  an  integrator  to  produce  a  passive  approximation  of 
a  zero  order  hold  circuit.    The  Pade1  polynomials  are  a  promising  method 
for  finding  a  rational  approximation  for  any  irrational  function  and  work 
can  be  done  in  the  implementation  of  such  functions.    Another  area 
where  work  can  be  done  is  to  consider  phase  relationships  in  the 
frequency  domain  analysis  of  the  error.     Only  the  amplitude  relationship 

(integral  square  error)  was  considered  in  this  thesis.   It  would  be 
interesting  to  see  what  additional  information  could  be  obtained  from 
such  an  analysis. 
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APPENDIX  1(c) 

-sT 
PADE'  APPROXIMATION  OF  1-e 

(Note:  all  listed  complex  roots  have  complex  conjugates  which  exe 
not  listed  but  are  represented  by  the  open  parenthesis  <> ) 
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G      =  X(x+3.6o^,u.577)(x+l.Uoo+.i3.700)(x-l.5QtH-,lU.306)()()n 

70  Tx+2.759)  (x-lTUo^i-jU.225)  (X+1.1U7+J3.12U)  (X+2.380+JI.629TDDI ) 

G       =  XjX^gjilgiU^^ 

71  (xT3TB39TTx-.509+(35.^95)  (x+2. 1*48+33. 950)  (x+3.Ma+j2.oi2)()(70 

G      =  X  (xr_._2i£+j^^^  . 6oU+ti7--338)  ( ) ,0-(  L....  v 

72  Tx?+7878)  (x-t-.UJ45-fj6.837)  (X+U. 1+70+32.325)  (x+3.l^SfjU76i6TDDD 


G    '=  X(X-.032+.l6.2U0)  (X+2.862+.19.61U)  (X-kUq170-i.15.H38)  ()()() 

73  Tx+T. 90il  (x+l .  iu8+ j  7 .  809 )  (X+5 .  ^Q6+t  2 .  598 )  (X+5 . lMi+j  5 . I9r+TDDD 

G       =  X (X+.012+ j 6 . 287 )  (X+2 . 500+.111 . 596)  (X+lU . 988+.18 . 625 1 ( )  ( )  0 

7h  TX+5T56B)  (x+2 . 206+J8 . 732 )  (x+U . 928+ j h . 726")  (x+7 . 532+j U .  5 67TDDT) 

G       =  XjX+iOOl+J^gtellXi^^ 

75  (x+7.928+j3.07l)(x+7.506)(x+3.338+j9.^50)(X+6.iU2+o6.i8U)()()() 

G       =  X  (X+ .  OOOU+j  6 .  283 )  (X+ .  508+.1 13 .  19M  (X+20 .  l+92+j  28 . 92MQO  0 

76  Cx^7?57)Tx1^^^  623TOOO 

G       =  X(X+tl6.283)(X+l1l2.958)(X+,136.n9)()()() 


=    X(X;M6.283)(X+lll2.958Hx+,136.n9)()()( 
CX+9.9*+3)(X+9  •  5l6^3TVf8)  (X+5 .371+jlo . 


77       Tx+9T^37fx-i^5l5+J3Tw)  (X+5O7l+jl0.8^l)  (x+8.ii!0+j7.03^)  00  0 
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